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Abstract 

Recently we proposed a new approach to the testing of dark energy models based on the ob- 
servational data. In that work we focused particularly on quintessence models for demonstration 
and invoked a widely used parametrization of the dark energy equation of state. In this paper we 
take the more recent SN la, CMB and BAO data, invoke the same parametrization, and apply this 
method of consistency test to five categories of dark energy models, including the ACDM model, 
the generalized Chaplygin gas, and three quintessence models: exponential, power-law and inverse- 
exponential potentials. We find that the exponential potential of quintessence is ruled out at the 
95.4% confidence level, while the other four models are consistent with data. This consistency test 
can be efficiently performed since for all models it requires the constraint of only a single parameter 
space that by choice can be easily accessed. 
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I. INTRODUCTION 



Compelling evidences from Type la supernovae (SN la) and other cosmological observa- 
tions show that the expansion of the universe is undergoing an accelerating stage at late 
times (see Ref. [l for reviews). Within the framework of general relativity and assuming 
homogeneity and isotropy, this indicates that there should exist an energy source, termed 
dark energy, which provides a significant negative pressure to cause this acceleration. Thus 
far the nature of dark energy remains unresolved and is generally regarded as one of the 
most tantalizing problems in cosmology. While a positive cosmological constant remains 
the simplest realization of dark energy, current observational data have not ruled out the 
possibility of a time-evolving dark energy, As one awaits more information from the future 
observations, the constraining power of the next-generation observations and new analysis 
methods are being pursued. j^-[lo| 

Many dark energy models have been proposed and studied. For the cosmological constant, 
its value has been constrained by observations (see Re f. I 111 , for example). The quintessence 
model, which invokes a time- varying scalar field. [f^-jlJl generally allows its energy den- 



sity and equation of state to evolve wit 
with different potential forms (see, Refs. 



Studies of the classification 



19|, 



;ime. There are various quintessence models 
for example) that have been proposed. 
20| and the general dynamical behavior 2]| of quintessence 



have been carried out. There have been also works on the reconstruction of quintessence 
potentials 21]- 28|,Q] and the investigations on how future observational data can constrain 



individual models of quintessence 



The generalized Chaplygin gas (see Re 



references therein) has been proposed to either unify dark matter and dark energy 

dark energy jsil, 



simply play the role o 
has been obtained. 



32] 



29 and 



30j or to 



32l ] . The constraint of the generalized Chaplygin gas 
3j] While these works have helped us study the possible nature of 
dark energy and obtain constraints of the parameters of an individual dark energy model, it 
should be desirable to explore possible means to determine whether a particular dark energy 
model can be ruled out by the observational data. 

Recently, we introduced a new approach to testing the consistency between a dark energy 



model and the observational data. 



8| We tested the exponential (see Ref. a and references 

n n 

therein) and the power-law potentials [1 71. Il8l ] of quintessence to demonstrate this method. 



which can be summarized as follows. For each dark energy model, we look for a charac- 



2 



teristic, Q{z), which in general can vary with the redshift but is equivalent to a constant 
parameter within the domain of the model. We further define the measure of consistency, 
A4{z), as the derivative of Q{z) with respect to the redshift z. The observational data 
should allow a null value for Ai{z) if the corresponding dark energy model is consistent with 
them. If, however, the A4.{z) = line lies outside certain confidence region, then that dark 
energy model is ruled out by the observational data at the corresponding confidence level. 
To obtain the constraint on the measure of consistency A4{z) from the observational data, 
a parametrization of the relevant physical quantity, such as the equation of state or the 
luminosity distance, is required. We have invoked a broadly used form of parametrization 



of the equation of state, 



35, M 



W{z) = Poe{z) / Poe{z) =Wo + Wa{l-a) =Wo + WaZ/{l + Z) , 



(1) 



where Pue{z) and Poe{z) are the pressure and the energy density of dark energy, respec- 
tively. The two parameters in Eq. ([T]) and the normalized matter density at present, Qm, 
define the parameter space, {wQ,Wa,^m), through which the information from the obser- 
vational data can be extracted. A recent work that is close in spirit to ours is that 



of Zunckel and Clarkson. 



10| who proposed consistency test of the cosmological constant 



via a direct parametrization of the luminosity distance. In our terminology, they use 
Q{z) = 1 — Poe{z) / Pc which is equivalent to the constant flm in the domain of the cosmo- 

I I 

logical constant, where pc is the critical density at present. Sahni et al.[9] also proposed null 
test of the cosmological constant via the diagnostic Om{z) = {H'^{z)/Hq — 1)/[(1 + z)^ — 1], 
which is equivalent to the constant in the domain of the cosmological constant, where 
H{z) is the Hubble expansion rate and Hq is the Hubble constant. 

In this paper, we take the more recent data set and apply our method of consistency test 
to five dark energy models, including the cosmological constant, the generalized Chaply- 

32| and three quintessence models: exponential, 



gin gas as the dark energy component. 



31 



power-law and inverse-exponential potentials 



17 



ISj. The data set we use includes a re- 



cently compiled "Constitution set" of SN la data,[37|-[43j the cosmic microwave background 
(CMB) measurement from the five-year Wilkinson Microwave Anisotropy Probe (WMAP) 
observation, [l])] and the baryon acoustic oscillation (BAO) measurement from the Sloan 
Digital Sky Survey (SDSS)^ and the 2dF Galaxy Redshift Survey (2dFGRS) 45|. 

A conventional way to determine how well a dark energy model can fit the observational 
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data is the model-based approach, in which one optimizes the parameters of each dark energy 
model based on the observational data and then statistically assesses the goodness of fit (see 
Ref. |3J, for example). In such approach one has to obtain the best fit for each set of param- 
eters specific to the particular dark energy model, which could be tedious. In particular, in 
order to optimize the parameters of a quintessence model one has to solve the field equation 
numerically for each point in the parameter space, which can be computationally intensive 
and time consuming, [^y] In contrast, in our approach we first constrain the parameters of the 
chosen parametrization through the observational data, and then test consistency of each 
dark energy model based on this set of parameters. This can be more efficient than the 
model-based approach when one deals with a large number of dark energy models. It is 
also more direct and therefore much faster to constrain the parameter space {wq, Wa, ftm) 
than to optimize the parameters of a quintessence model in the model-based approach. The 
potential downside of our method, however, would be that as long as one invokes a specific 
form of parametrization, one might have simultaneously imposed a prior, or bias, against 
certain dark energy models. This issue requires a separate investigation and we are currently 



pursuing that.[46|] We note that the two methods are different in spirit. The goodness of 
fit describes how well a model can fit the observations. The consistency test, on the other 
hand, examines whether the condition necessary for a model is excluded by the observations. 
With in mind the pros and cons mentioned above, we believe that the two methods, that is, 
the model-based and ours, should be complimentary to each other in the pursuit of revealing 
the nature of dark energy. 

II. CONSISTENCY TEST OF DARK ENERGY MODELS 
A. Formalism 

In this paper we consider a fiat Friedmann-Lemaitre-Robertson- Walker (FLRW) universe 
and assume that it is dominated by pressureless matter and dark energy in the present 
epoch. The Hubble expansion rate, H = a/a, is given by the Friedmann equations as 



H {z) = [pm{z) + Pde{z)\ 



Hi fi„(l + ^)3 + (l-fiJexp(^3^ [1 + «;(/)] 



(2) 
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where the dark energy density 

Pde{z) = pc{l - Qm) [1 + Yq^^ , (3) 

and 

ottG at 

For quintessence as a dark energy model, the quintessence field and the potential are related 
to the equation of state, the Hubble expansion rate, and the dark energy density as follows 

0,0. 

Viz) = [l-wiz)]p^^{z)/2. (6) 



We perform the consistency test of five dark energy models including the cosmological 
constant, the exponential potential, the power-law potential, the inverse-exponential poten- 
tial, and the generalized Chaplygin gas as the dark energy component. 

For the cosmological constant, the energy density pA is a constant. We define the char- 
acteristic Qa{z) as the dark energy density Pde{z), which in general would evolve with the 
redshift but is equivalent to the constant parameter p\ within the cosmological constant 
domain, 

Qa(2) = Pde(^) (7) 
= Pa for the cosmological constant. (8) 

In the same spirit, for the exponential potential, 

V;xp(0) = V^iexp[-0/Mi] , (9) 

we identify Mi as the characteristic constant parameter and accordingly define the charac- 
teristic Qexp{z), 

Q^^,{z) ^ -Viz) ' (;.) (10) 

= Ml for the exponential potential. (11) 

For the power-law potential, 

V^power(0) =m^->", (12) 
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we define the following characteristic corresponding to the index n, 



Qpower(-2') 



1 T/r ^ /^^^^ ^^"'^'^ 



n for the power-law potential. 



For the inverse-exponential potential, 



mverse-exp 



(0) = l/2exp[M2/(/)] , 



(13) 
(14) 

(15) 



the characteristic is defined as 



Q 



mverse-exp 



V{z) v# 



dV 



d'^V , , 1 



V{z) \d(t) 



dV 



M2 for the inverse-exponential potential. 



(16) 
(17) 



As the dark energy component, the generalized Chaplygin gas has an equation of state 
govern by 

p^z) ^ -A/ [p,,iz)r , (18) 
where a —1 and ^4 > 0. The corresponding characteristic is defined as 



Q 



Chaplygin 



w{z) dz \ dz 
— a for the generalized Chaplygin gas. 



(19) 
(20) 



We then define the measure of consistency M.i{z) as the derivative of the characteristic 
Qi{z) with respect to the redshift for each dark energy model, 

dQi 



MAz) 



dz 



for the corresponding dark energy model. 



(21) 
(22) 



where i denotes "A", "cxp", "power", "inverse-exp" and "Chaplygin", respectively. M.i{z) 
can in general evolve with the redshift but should be constant zero in the domain of the 
corresponding dark energy model. 



B. Observational data and constraint 



We use the combined data set from three types of observations including the SN la 
observation, the CMB measurement, and the BAO measurement. 
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We use the Constitution set of SN la data compiled by Hicken et al., 43|, 37|]- 42] which 
provides the information of the luminosity distance and the redshift. The luminosity 
distance-redshift relation is given by 

dz' 



dr.(z) 



l + z) 







H(z' 



We use the CMB shift parameter measured by the five-year WMAP observation, 

.1090.04 



R 



n H"^ 



1.710 ±0.019. 



We use the BAO measurement from the joint analysis of the SDSS and 2dFGRS data, 
which gives 



(23) 



(24) 



44 



45| 



Dy(0.35)/Dy(0.2) = 1.812 ± 0.060, 



where 



Dviz 



:i±^B 



Diiz 



Zbko 



Hiz, 



BAOj 



and Da{z) is the angular diameter distance, 

Da{z) 



l + z 



dz' 
Hiz' 



1/3 



(25) 
(26) 

(27) 



The constraint of the parameter space (wq, Wa, fi^) is obtained by fitting the three 
parameters to this combined data set. The best fit of the parameters are found to be 



Wo 



-0.89l°:}2, wa = -0.18 



+0.71 
0.745 



0.25 



+0.03 
-0.02- 



The two-dimensional constraint of {wo,Wa) is shown in Fig. [H 



C. Results of the consistency test 

For the consistency test of each dark energy model, we reconstruct Adi{z) via the con- 
straint of {wq, Wa, ^m), with the use of the equations in Sec. [land Sec. Ill A[ We perform the 
test in the redshift region < z < 1.55, where the influence of dark energy on the expansion 
of the universe is most significant. This region is covered by the current SN la observations, 
which is the most sensitive type of observations to probe the behavior of dark energy. If the 
Aii{z) = line lies outside certain confidence region, the corresponding dark energy model 
is ruled out at that confidence level. Adopting the constraint obtained in Sec. IIIB[ we find 
that the Aicx.p{z) = line lies outside the 95.4% confidence region while the null lines of the 
other four models lie inside the 68.3% confidence region. The results are shown in Fig. [2j 
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III. SUMMARY 



We have preformed consistency test of five dark energy models, including the cosmolog- 
ical constant, the generalized Chaplyngin gas, and three quintessence models: exponential, 
power-law, and inverse-exponential potentials. This test gives a simple signature if a dark 
energy model is ruled out by the observational data. It can be done efficiently via the 
constraint of a single set of parameters deduced from the observational data, and can test 
quintessence models without solving the field equation. 

Through our approach and invoking the broadly used parametrization of the equation of 
state, the exponential potential is found to be ruled out at the 95.4% confidence level based 
on the current observational data. The other four dark energy models remain consistent 
with the current observations down to the 68.3% confidence level. It is worth noticing that 
in our previous work the power-law potential was ruled out at the 68.3% confidence level 
based on a different data set [8]. One noticeable change in the new data set is that the 
cosmological constant is contained in the 68.3% confidence region, which was not so in the 
previous one. Whether our method can discriminate between the power-law potential and 
the cosmological constant can in principle be studied with the Monte Carlo test 46|]. The 
fiat cosmological constant model was examined by Davis et al ^]. via the model-based 
approach, in which they used the same data set except the three-year WMAP data instead. 
They found the goodness of fit for the fiat cosmological constant model to be 43.7% while 
we find the model at least consistent with the data at the 68.3% confidence level. The two 
results are not in confiict with each other. 

Our method of the consistency test can in principle be applied not only to other dark 
energy models but also to other models explaining the acceleration of the expansion, as 
long as one looks for the characteristic Q{z) corresponding to a constant parameter of each 
model. One can also choose a different parametrization for better discriminating power 
between the models in regard. The discriminating power of the method with different forms 
of parametrization and the possible bias imposed by the chosen parametrization should be 



studied via the Monte Carlo test 



46|. 
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FIG. 1: The joint two-dimensional constraint of {wo,Wa) based on a combined data set including 
the Constitution set of SN la data, the CMB measurement from the five-year WMAP, and the 
BAO measurement from the SDSS and 2dFGRS. The dark and light gray areas correspond to the 
68.3% and 95.4% confidence regions, respectively. 
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FIG. 2: The measure of consistency of the five dark energy models. The dark and light gray areas 
correspond to the 68.3% and 95.4% confidence regions, respectively. The Mexp{z) = line lies 
outside the 95.4% confidence region for 0.95 < z < 1.55. The null lines of the measure for the 
other four models lie inside the 68.3% confidence regions for < z < 1.55. This indicates that the 
exponential potential is ruled out at the 95.4% confidence level while the other four dark energy 
models are still consistent with the current observational constraints down to the 68.3% confidence 
level. 
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